International Journal of Mechanical and Production 
Engineering Research and Development (IJMPERD) 
ISSN(P): 2249-6890; ISSN(E): 2249-8001 
Vol. 9, Issue 4, Aug 2019,307-316 
© TJPRC Pvt. Ltd. 



TRAN5 

STELLAR 

• Journal Publications • Research Consultancy 


THERMAL ANALYSIS FORINTERNAL COMBUSTION 
CYLINDER BLOCK 

MUBARAK F. AL-DHAFIRI & HASAN A. ALZENKI 

The Public Authority for Applied Education and Training 
High Institute ofEnergy, Kuwait 

ABSTRACT 

This paper presents a simulation and analysis for the thermal stresses of the internal combustion engine and 
analyzed its effects on cylinders and engine performance. The main aim of this paper is to make a simulation and 
analysis of the thermal stresses affect internal combustion engine or thermal fatigue. A simulation using Solidworks is 
used here to analyze the thermal stresses of ICE problem. It is found that the thermal stresses of cast carbon are of 
smaller value than that ofthe alloy steel. Thermal displacements and strains of alloy steel blocks are smaller than that of 
cast carbon. Cooling contributes in reducing thermal stresses, strains and displacements on both different blocks. 
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INTRODU CTION 

The thermal fatigue of internal combustion engines is the principle on which this research is based on the 
used engine in experimental tests was a Dacia engine. If the negative effects in the case of any damage to the gears 
that constitute internal combustion engines, in order to mitigate the destruction caused by cracking, some 
comprehensive, complex, theoretical and empirical research on thermal fatigue, using theoretical and experimental 
approaches should be conducted. The study was conducted in three phases: first studying the thermal component of 
the partial mechanisms, to study and conduct an experimental determination of the changing fields of the surface 
temperature of the components of the device, and to perform an experiments on heat life samples exposed to 
different forms of thermal stress. 

This study aimed to is to solve the complex problems of internal combustion engines (ICE) by controlling 
a new, original, theoretical and empirical knowledge about the heat exhaustion occurring in machines, which 
allows the transfer of heat from the combustion chamber to the cooling fluids of motors installed on automobiles 
[ 1 ]. 

It is known that cracks resulted from thermal fatigue appear on the surface and inside the outer layers of 
heat transfer parts. The specific thermal fissures develop gradually due to periodic changes in temperature. 
These cracks, on the upper layers and on the surface of the machine parts (although operating conditions are most 
appropriate), limit their use in situations that have adverse effects on the energy and economic indicators of the 
engine [2], [3]. 
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Figure 1: Experimental Bench to Measure the Temperature Values of ICE 

There are many research discussed this issue; in this study, the researcher tried Amit V. Paratwar, D.B. Hulwan 
(2013) heat transfer and process flow analysis. The cover of the current average diesel engine cooling turbine refrigeration 
works of 6 cylinders, and examine the factors that affect cooling performance. To optimize the parameters listed CFD 
analysis of steady state verification and validated experimental results. 

Martin G. and S. Olle (2015) tried to develop a model of diesel piston thermal analysis was used in the engines of 
Volvo Cars. Enhances thermal analysis model for identifying potential piston thermal loads full of drive and can be used as 
condition limits when CFD simulation of combustion. The study might be a help in the design of the piston. 

Shobam. And Shekhar (2016) performed thermal analysis of engine block with fins. Parametric model was 
developed for the engine block fins in a 3D program Solidworks is made and thermal analysis on the fins and mass to 
determine the temperature change in the interim situation in stable condition with and without looking over time. 
The analysis is performed by using different materials. In this study two models were created in the program and the 
modified design was analyzed the form itself, and a comparison of the two models according to the geometry and the 
materials were analyzed. New materials were selected such as aluminum alloy 1050 for this analysis. 

Vinay A. and Himanshu A.(2016)a three dimensional solid model of piston including piston pin is designed with 
the help of Solidworks software was constructed. The temperature, total heat flux are calculated to investigate factor of 
safety and life of the piston assembly using ANSYS workbench software. Aluminum alloy was used as piston material. 
The transient stress analysis results also help to improve component design at the early stage and also help in reducing time 
required to manufacture the piston component and its cost. 

Vinay A. and Himanshu A.(2016) constructed a 3D solid model design of the piston, including the piston pin with 
the help of Solidworks software. The temperature is calculated, the total fiow Heat and Safety factor and piston set age 
using ANSYS Workbench. Aluminum alloy is used as a piston. Stress test results also help transit to improve the design of 
the components in the early stage and also help to reduce the time required to manufacture the piston and component cost. 

MODELLING AND SIMULATION 

The modelling of the ICE is performed using the Solidworks software application. Figure 2 shows the cylinders of 
the ICE and engine block and its meshing. This engine consists of four cylinders. All specifications of the model and 
modelling parameters. 
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Figure 2: Cylinder Block of the ICE and its Meshing 

This simulation will show engine thermal stress inside cylinders in every stroke, the simulation will be in two 
different ways, the first one is showing thermal stresses without engine cooling, and the second one will be with engine 
cooling, water convection data as in table 1 below. 

Engine firing order starts from cylinder number 1,3,4,2 ( Red, Brown, Blue, Light Blue) so the thermal stress on 
each cylinder is shown on the Table 1 below, the cooling process for the engine is applied to show cooling efficiency. 


Table 1: Thermal Loads 


Load Name 

Load Image 

Load Details 

Red 

Cylinder 





Entities: 

Temperature: 

Power Stroke 

500C 

Blue 

Cylinder 


■)' * \ 



Entities: 

Temperature: 

Intake Stroke 

37C 

Brown 

Cylinder 


- m 

X. rJ 

;( M> : ; 



Entities: 

Temperature: 

exhaust stroke 

200 C 

Light Blue 
Cylinder 


jffif 



Entities: 

Temperature: 

compression stroke 

110C 

Engine Head 





Entities: 

Convection 

Coefficient: 

Engine Head cooling 
from water 

1000 W/(m A 2-K) 
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RESULTS AND DISCUSSIONS 


Figure 3 shows the temperature distribution with cooling of cast iron engine block. The red areas have more 
temperatures than that of blue areas. 


Model name SW3dPS-4 cydnder crank assemblyl 
Study nemet Ergne coolrig 
Plot type: Thermal Thernml! 





Temp (Cetseis) 

5 000e»002 
L 4,614e+002 


I 


. 3.842e*002 

. 3 456ee002 
. 3.0710*002 


2885e*C02 
2.299e*002 
18136*002 
1 527e*002 


1.1418*002 


I 


75516*001 

36926*001 


Figure 3: Temperature Distribution of Cast Carbon ICE Block 


Figure 4 shows the Cast carbon ICE block thermal stress-Von-Mises stresses. 


Model nemr SW3dPS-4 cydnder crank assemfc(y2 

Studyname static 

Plot type Statlc nodal stress Stress! 


7l 



I 


ases (N/ln A 2) 
99t ,662,200 0 
909062.7840 
826.463,3600 
743063,9380 
681064,5120 
578066,0880 
496 065,684 0 
413,4682400 
330068,8160 
248267,3920 
165067,9880 
83068,5520 


♦ VieWstrengttv 248,168000.0 


Figure 4: Cast Carbon ICE Thermal Stresses Distributions 


Figure 5 shows the displacement of the cast Carbon ICE Block. The red areas have more displacements and 
thermal strains. Also figure 6 shows the thermal strains of the same block. 


Modd name: SW3dPS-4 cyclnder crar* essembly2 
Study name. stabc 

Ptot type Siatic cisptacemert Displacementt 



Figure 5: Displacements of Cast Carbon Engine Block 
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Figure 6 shows the thermal strains distribution of the cast carbon ICE block. 


Modd name: SW3dPS-4 cyclinder crank assemb<y2 

Study name: static 

Plot type: Static stran Strainl 



Figure 6: Thermal Strains Distribution of the Cast Carbon ICE Block 

Figure 7 shows the cast carbon ICE cylinder block temperature distribution without cooling. 


Model name SW3dPS-4 cycSnder crank assembly2 
Study name [Englno no coorngj 
Plot type: Thermal Thermall 
Time stepr 1 





Temp (Celsius) 

1 5 0006+002 
4.614e+002 
4 226e+002 
. 3 843e+002 
. 3 457e+002 
. 30716+002 


. 22996+002 
. 10136+002 
1 528e+002 
. 1 1428+002 
r 7 5S8e+001 
l 3700e+001 


Figure 7: The Cast Carbon ICE Cylinder Block Temp. Distribution without Cooling 


Figures 8 to 10 show the Cast carbon block: thermal stress, strain and displacement. 


Model name SW3dPS-4 cydnder crank assembiy2 

Study name: static no cooling 

Ptottype Static displacemert Displacementl 



I 


9843e-001 
89486-001 
. 8053e-001 
. 71586-001 
. 6 264e-001 
. S569e-001 
. 4 4746-001 
3 5796-001 
2 684e-O01 
. 1 790e-001 
t 89486-002 
1 1 0006-030 


Figure 8: The Displacement of the Cast Carbon Block without Cooling 
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Modd name SW3dPS-4 cyctnder erank assembty2 
Studyname static no coobng 
Plot type Static straln Strani 



ESTRN 



. 4 485e-003 

. 3 968e-003 
. 3 491e-003 
. 2S93e-003 
. Z496e-003 


. 1 501e-003 
. 1 004e-003 

[ 5fl67e-004 
9.410e-006 


Figure 9: The Strain Distribution of the Cast Carbon Block without Cooling 


Model name: SW3dPS-4 cyclinder crank assembly2 

Study name: static no cooling 

Plot type: Static nodal stress Stressl 



von Mises (NAn A 2) 

P 1,888,269.312.0 
1,731,044,096.0 

_ 1,573,818,752.0 

. 1,416,593,536.0 
. 1,259,368,320.0 
. 1,102,143,104.0 
. 944,917,824.0 
. 787,692,544.0 
. 630,467,328.0 
. 473,242,112.0 
. 316,016,864.0 
bd 158,791,616.0 
■ 1,566,364.0 
-► Yield strength: 248,168,000.0 


Figure 10: The Thermal Stresses Distribution of the Cast Carbon Block without Cooling 

Figure 11 shows the alloy steel ICE block temperature distribution with cooling. 


Model name: SW3dPS-4 cyclinder crank assembly2 
Study name: SS COOLING 
Plot type: Thermal Thermall 
Time step: 1 



Temp (Celsius) 

I 5.000e+002 
4.614e+002 
4.228e+002 
. 3.842e+002 
. 3.457e+002 
. 3.071 e+002 
. 2.685e+002 
. 2.299e+002 
. 1.913e+002 
. 1.527e+002 
. 1.141e+002 
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■ 3.697e+001 


Figure 11: Alloy Steel ICE Block Temperature Distribution with Cooling 
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Figure 12 to 14 show the alloy steel thermal stress, strain and displacement of the ICE with cooling. 


Model name: SW3dPS-4 cyclinder crank assembly2 
Study name: static SS with cooling 
Plot type: Static nodal stress Stressl 



von Mises (Nftn A 2) 


1 1,283,867,904.0 
1,176,969,472.0 
_ 1,070,071,232.0 
. 963,172,928.0 
. 856,274,560.0 
. 749,376,192.0 
^ . 642,477,888.0 

. 535,579,520.0 
. 428,681,184.0 
. 321,782,848.0 
. 214,884,480.0 
bJ 107,986,136.0 
H. 1,087,794.3 


-► Yield strength: 620,422,000.0 


Figure 12: Alloy Steel Thermal Stress of the ICE with Cooling 


Model name: SW3dPS-4 cyclinder crank assembly2 
Study name: static SS with cooling 
Plot type: Static strain Strainl 



ESTRN 


I 3.507e-003 
3.216e-003 
2.924e-003 
. 2.632e-003 
. 2.340e-003 
. 2.048e-003 
. 1.756e-003 


. 1.173e-003 
. 8.808e-004 
. 5.890e-004 
jl 2.971 e-004 
I 5.292e-006 


Figure 13: Alloy Steel Thermal Strain of the ICE with Cooling 


Model name: SW3dPS-4 cyclinder crank assembly2 
Study name: static SS with cooling 
Plottype: Static displacement Displacementl 


URES (mm) 

| 7.366e-001 
6.752e-001 
— 6.138e-001 
. 5.525e-001 
. 4.911e-001 
. 4.297e-001 
. 3.683e-001 
. 3.069e-001 


. 2.455e-001 
. 1 842e-001 
. 1.228e-001 
J 6.138e-002 
L 1 000e-030 


Figure 14: Alloy Steel Thermal Displacement of the ICE with Cooling 
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Figure 15 shows the alloy steel ICE block temperature distribution without cooling. 


Model name: SW3dPS-4 cyclinder crank assembly2 
Study name: SS NO COOUNG 
Plottype: Thermal Thermall 
Time step: 1 





Temp (Celsius) 

I 5.000e+002 
4.614e+002 
4.228e+002 
. 3.843e+002 


. 3.457e+002 
. 3.071 e+002 


. 2.685e+002 
. 2.299e+002 
. 1.913e+002 
. 1.528e+002 
. 1.142e+002 
l 7.558e+001 
L 3.700e+001 


Figurel5: Alloy Steel ICE Block Temperature Distribution without Cooling 

Figure 16 to 17 show the alloy steel thermal stress, strain and displacement of the ICE block without cooling. 


Model name: SW3dPS-4 cyclinder crank assembly2 
Study name: Static SS no cooling 
Plot type: Static strain Strainl 



ESTRN 

| 6.279e-003 
5.757e-003 
. 5.234e-003 
. 4.712e-003 
. 4.189e-003 
. 3.667e-003 
. 3,144e-003 
. 2.622e-003 
. 2.100e-003 
. 1.577e-003 
. 1.055e-003 

y ‘ J 5.322e-004 
9.737e-006 


Figure 16: Alloy Steel Thermal Strain of the ICE Block without Cooling 


Model naira SW3dPS-4 cyctnder crar* sswetv2 
Shrdy naire Stalic SS no codrg 
PM type Stitc nodal streas Stressl 




I 


Ms«(Niln*2) 
2,155.4935860 
1.9rejJ11,520.0 
1.796,529024 0 
. 1817jO«5280 
. 1,4375640320 
. 1,258fl81536.0 
10785990400 
899.1165080 
719JB34.1120 
. 540,1510160 


1 181.1865860 
1.7045114 

+ vield arengtn 620,4220000 



URESdnm) 



. 6 723*001 
. 77S4e-001 

. 6 784*001 
. 5815*001 
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(a) (b) 

Figure 17: Alloy Steel Thermal Stress and Displacement of the ICE Block without Cooling 


Impact Factor (JCC): 7.6197 


SCOPUS Indexcd Journal 


NAAS Rating: 3.11 




Thermal Analysis for Internal Combustion Cylinder Block 


315 


CONCLUSIONS 

Studying the thermal stresses and temperature distributions and also the thermal strains and displacement of the 
ICE block with different materials is very important issue especially using qualified finite and simulation and modelling 
software like Solidworks gives the opportunity to determine which of the materials studied is the best. It is found that the 
thermal stresses of cast carbon are of smaller value than that of the alloy steel. Thermal displacements and strains of alloy 
steel blocks are smaller than that of cast carbon. Cooling contributes in reducing thermal stresses, strains and displacements 
on both different blocks. 
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